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Abstract. The aim of this Note is to prove by a perturbation method the 
existence of solutions of the coupled Einstein-Dirac-Maxwell equations for a 
static, spherically symmetric system of two fermions in a singlet spinor state 
and with the electromagnetic coupling constant ( — ) 2 < 1. We show that the 
nondegenerate solution of Choquard's equation generates a branch of solutions 
of the Einstein-Dirac-Maxwell equations. 

Une methode de perturbation pour les solutions localisees des equa- 
tions d'Einstein-Dirac-Maxwell. 

Resume. Le but de cette Note est de demontrer par une methode de pertur- 
bation l'existence de solutions des equations d'Einstein-Dirac-Maxwell pour 
un systeme statique, a symetrie spherique de deux fermions dans un etat de 
singulet et avec une constante de couplage electromagnetique ( — ) < 1. On 
montre que la solution non degeneree de l'equation de Choquard genere une 
branche de solutions des equations d'Einstein-Dirac-Maxwell. 



Version franqaise abregee 

Dans un papier recent [5], par une methode de perturbation, on a montre de 
maniere rigoureuse l'existence de solutions des equations d'Einstein-Dirac pour un 
systeme statique, a symetrie spherique de deux fermions dans un etat de singulet. 
Dans cette Note, on generalise ce resultat aux equations d'Einstein-Dirac-Maxwell 
et on montre, dans le cas particulier d'un couplage electromagnetique faible, l'exis- 
tence des solutions obtenues numeriquement par F. Finster, J. Smoller et ST. Yau 
dans 0. 

Plus precisement, en utilisant l'idee introduite par Ounaies pour une classe 
d'equations de Dirac non lineaires (voir [1]) et adaptee dans [5] aux equations 
d'Einstein-Dirac, on obtient le theoreme suivant. 

Theoreme 0.1. Soient e,m,u> tels que e 2 — m 2 <0,0<uj<met supposons 
m — lo assez petit ; alors il existe une solution non triviale de {7]|5[). 

Dans cette Note, on decrit la methode utilisee pour demontrer ce theoreme. 
Premierement, par un changement d'echelle, on transforme les equations d'Einstein- 
Dirac-Maxwell ([6][9]) en un systeme perturbe qui s'ecrit sous la forme (fTTjl . On choisit 
e = m — uj comme parametre de perturbation. 

Deuxiement, on remarque que, pour e — et (^) < 1, ce systeme est equivalent 



au systeme Ifl2|) ou l'equation pour la variable ip est l'equation de Choquard. II est 
bien connu que l'equation de Choquard a une solution radiale positive. De plus, dans 
l'espace des fonctions radiales, cette solution est non degeneree, dans le sens ou le 
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noyau de la linearisation de l'equation contient seulement la fonction identiquement 
nulle. On appelle <po la solution du systeme ([12]) . 

Ensuite, on observe que le systeme perturbe s'ecrit sous la forme D{e, ip, x, t, Q = 
avec D un operateur non lineaire de classe C , pour un bon choix d'espaces 
fonctionnels. On prouve que cet operateur satisfait les hypotheses du theoreme des 
fonctions implicites. En particulier, on montre que la linearisation de Poperateur 
D par rapport a (<p,X> T > en (0><^o); D VyX:T ^(0, (j> ), est une injection, grace a 
la non-degenerescence de la solution de l'equation de Choquard, et s'ecrit comme 
somme d'un isomorphisme et d'un operateur compact ; done D v ^ XiT ^(0, 4>q) est un 
isomorphisme. En appliquant le theoreme des fonctions implicites, on deduit que, 
pour e assez petit et e 2 — m 2 < 0, le systeme (fTTj) a une solution. 

En conclusion, pour e 2 — m 2 < 0, < lj < m et m — u> assez petit, les equations 
d'Einstein-Dirac-Maxwell possedent une solution non triviale. 



1. Introduction 

In a recent paper [5], using a perturbation method, we proved rigorously the ex- 
istence of solutions of the coupled Einstein-Dirac equations for a static, spherically 
symmetric system of two fermions in a singlet spinor state. In this Note, we extend 
our result to the Einstein-Dirac-Maxwell equations and we prove, in the particular 
case of a weak electromagnetic coupling, the existence of the solutions obtained 
numerically by F. Finster, J. Smoller and ST. Yau in [7]. 

The general Einstein-Dirac-Maxwell equations for a system of n Dirac particles 
take the form 



(G - m)ip a = 0, R) - lR8) = -8nTj, W k F^ k = Aire £ ^ a G^, 



where are the Dirac matrices, G denote the Dirac operator, ip a are the wave 
functions of fermions of mass m and charge e, Fjk is the electromagnetic field tensor 
and, finally, Tj is the sum of the energy-momentum tensor of the Dirac particle and 
the Maxwell stress-energy tensor. 

In [7], the metric, in polar coordinates (i, r, z9, i^), is given by 

ds 2 = \dt 2 -— dr 2 - r 2 dti 2 - r 2 sin 2 d dip 2 



with A = A(r), T — T(r) positive functions; moreover, using the ansatz from [6], 
Finster, Smoller and Yau describe the Dirac spinors with two real radial functions 
<I>i(r), <I>2(r) and they assume that the electromagnetic potential has the form 
A = (—V, 0), with V the Coulomb potential. 
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In this case the Einstein-Dirac-Maxwell equations can be written as 
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with the normalization condition J °° |3| 2 ^7= dr = 

In order that the metric be asymptotically Minkowskian and the solutions have 
finite (ADM) mass, Finster, Smoller and Yau assume 

lim T(r) = 1 

r — >oo 

and 

T 

lim -(1 - AM) < 00. 

Finally, they also require that the electromagnetic potential vanishes at infinity. 

In this Note, using the idea introduced by Ounaies for a class of nonlinear Dirac 
equations (see [T]) and adapted in [5] to the Einstein-Dirac equations, we obtain 
the following result. 

Theorem 1.1. Given e,m,oj such that e 2 — m 2 < 0, < w < m and m — u is 
sufficiently small, there exists a non trivial solution of (Tj^j). 

2. Perturbation method for the Einstein-Dirac-Maxwell equations 

First of all, we observe that, writing T(r) = 1 + t(r) and integrating the equation 
(JSj) , the Einstein-Dirac-Maxwell equations become 

VI*; = i$ 1 -(( w -eV)(l + t)+m)$ 2 (6) 
r 

VA& 2 = ((w - eV)(l + 1) - m)$i - -<S> 2 (7) 

r 

2rAt' = (A - 1)(1 + t) - 16tt(w - eV)(l + tf ($ 2 + 

+327ri(l + t) 2 $i$ 2 + 167rm(l + i) 2 ($ 2 - $ 2 ) 

+r 2 A(l+t) 3 (^') 2 (8) 
^(l+t)V = -^jT(*» + *»)£+^d.. (9) 

where A(r) = 1 + a(r) and 

1 /"" 
o(r) = — exp(-F(r)) / [16tv(w - eV)(l + i) 2 ($ 2 + $|) 
r Jo 

+s 2 (l + i) 2 (y') 2 l exp(F(s))ds (10) 
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with F(r) = J r S(l + tf {V'f ds. 

After that, we introduce the new variable (<p, x, T , C) such that 

$ x (r) = e 1 /2 (/ ,( e i/2 r ) 5 $2 ( r) = £X ( £ V2 r ), t(r) = er^r), V{r) = eC(e 1/2 r) 

where $i,$2j^, V satisfy (JSHSJ) and £ = m — u>. Using the explicit expression of 
a(r), given in (fTO)) . we write 

a($i, $ 2 , i, U) = ea(e, ip, x, r, C) 

with a(0, i/3, x, r, () = — I^l™ p 2 ds . It is now clear that if $1, <j> 2 , t, V satisfy 
([6]l9l, then x, r, £ satisfy the system 

' (1 + ea(e, <p, X , t, C)) 1/2 &<p -$<p + 2m X + K x (e, x, r, = 
(1 + ea(e, p, x, r, 0) 1,2 ^X + mcpr + ep( + K 2 (e, cp, x, t,Q = 
(1 + ea(e, <p, X , r, C))£t - °(^*^0 + K 3 (e, <p, X) r, Q = 

t (1 + ea(e, <p, X , t, C)) 1/2 (l + er)££ + ^ f ds + K 4 (e, <p, X , t,Q = 

(11) 

where Ki (0, ip, x, r, C) = K 2 (0, <p, x, r, C) = ^3 (0, </J, X, T , C) = #4 (0, p, X, r, Q = 
0. 

Then, for e = 0, (jXTJ) becomes 

' + 2my> + 167r(e 2 - m 2 )m (/ °° ds) ip = 

X(r) = ^(^-^) 

r(r) = 87rmJ oo m J ( 2 rs) ds 
^ CM =87re/ °° m& f {r s) ds 

We remark that if e 2 — m 2 < the first equation of the system (fT2)) is the Choquard 
equation 

-Au + 2mu-4(m 2 - e 2 )m (/ R3 cfy) u = in i? 1 (R 3 ) (13) 

with u(a;) = , It is well known that Choquard's equation lfl3|) has a unique 
radial, positive solution w with / |w | 2 = N for some N > given. Furthermore, 
uo is infinitely differentiable, goes to zero at infinity and is a radial nondegenerate 
solution; by this we mean that the linearization of lfl3|) around uq has a trivial 
nullspace in L 2 (R 3 ) (see [5], 0], [2] for more details). 
Let 0o = (po, Xo, to, Co) be the ground state solution of (J2|) . 

The main idea is that the solutions of ifTTj) are the zeros of a C 1 operator D : 
R x X v x X x x X T x X{ ^> Y<p xY x xY T xY$. So, to obtain a solution of (|TT|) from 
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t>o, we define the operators 



Li(e,V,X.r,C) = (1 + ea{e, <p, X , r, 0) 1/2 -~r<P - + 2 



1 d ip x 



+ -K 1 (s,tp,x,T, C) 



L>2{e, <P,X, T , C) 
L3(£,<P,X,t, C) 

and 



dr ' 



(1 + ea(e, ip, x, r, 0) 1/2 ~X + 4 + " " + ^ 

r ar r r r r 

+-K 2 (s,tp,x,T, C) 
r 

(1 + ea[e, ip, X , t, Qj^t ^ 1- K 3 (e, f, X, t, Q 



, , d ^ 87re 

(1 + ea(e, ip, X , r, C)) 1/2 (l + er) — £ + — 

ar r z 

+K 4 (e, ip, x, t, C) 



ip ds 



D(e, ip, x, t, C) = (Li(e, ip, x, r, (), L 2 (e, p, x, r, Q, L 3 (e, ip, x, t, Q,L A (e, ip, x, r, C)) , 
with X v , X x , X T , Y ip ,Y x ,Y T denned as in [5] and 



X r 



C:(0,oo) 



lim CM =0,-f-C€ ^((0, oo), dr) ni 2 ((0,oo), rdr) \ 

r^oo Ctr J 



Y c = L 1 ({Q,oo),dr)nL 2 ({0,oo),rdr) 

with their natural norms. 

Next, we linearize the operator D on (<p, x, T, C) around (0, (j>o): 



( 



D v , x , T x(0,4>o)(h,k,l,z) = 



r dr r 



1 d 



ar r r r r 



■2-1 

dr 1 



I 



dr 

\ 

-m^T + e^Co 
^ Jo 'fiohds 
V ^J^ohds / 

We observe that, thanks to the nondegeneracy of the solution of Choquard's equa- 
tion, D ViXtT ^(0, 4>q) is a one-to-one operator. Moreover, it can be written as a sum 
of an isomorphism and a compact operator. It is thus an isomorphism. Finally, 
the application of the implicit function theorem yields the following result, which 
is equivalent to Theorem ll.il 

Theorem 2.1. Suppose e 2 — m 2 < and let 4>o be the ground state solution of 
(CUP, then there exists 5 > and a function 77 <E C((0, 5), X v x X x x X T x X^) such 
that 77(0) = 4> a and D(s, 77(e)) = for < e < 5. 
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